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THE REMOTE SENSING OF SATURN'S RINGS 


I: THE MAGNETIC ALINEMENT OF THE RING PARTICLES 

Lawrence C . Evans 
Ames Research Center 

SUMMARY 


Because of the potential implications for the optical properties of 
Saturn's rings, the orientation of nonspherical ring particles at equilibrium 
is investigated with respect to four stochastic influences: interactions with 

the interplanetary medium, interactions with the expected magnetic field of 
Saturn, thermal fluctuations due to the internal temperature of the ring par- 
ticles, and collisions between ring particles. The solution of the homoge- 
neous Fokker-Planck equation for nearly spherical spheroids is presented and 
investigated in general. Values of the pertinent physical parameters in the 
vicinity of Saturn are estimated, and the implications for the alinement of 
the ring particles are investigated. It is concluded that for some alinement 
mechanisms, small ring particles (size < 10 cm) can be expected to be almost 
completely alined (a < 1°) . This alinement results in each particle spinning 
around its shortest body axis with this axis parallel to the magnetic field 
direction (perpendicular to the ring plane) . 

INTRODUCTION 


This is the first in a series of studies dealing with the potential sci- 
entific data obtainable by remote sensing of Saturn's rings from a spacecraft. 
This study deals with the configuration of nonspherical ring particles that 
would be expected in the presence of a Saturnian magnetic field. The second 
study (Part II) 1 is an analysis of the geometrical considerations of such 
remote sensing and the interaction of these considerations with such items as 
photometer thresholds, phase angle coverage, and areal and radial resolution. 
Subsequent studies (e.g.. Part III) 1 will deal with the implications of the 
expected particle configuration for remote sensing data. 

This series of studies has been prompted by the growing interest in fly- 
by and orbiter missions to Saturn. Although the major thrust of the explo- 
ration of the outer planets is currently directed toward Jupiter, Saturn is 
the major secondary objective. One of the reasons for the sustained scien- 
tific interest in Saturn is undoubtedly the planet's rings. This interest 
is concentrated in two main areas: the features of the rings themselves, 

and the cosmological significance of their existence. The rings of Saturn 
are one of the few unique objects in the solar system; as such, the intrinsic 
interest in them is certainly not surprising. Their very existence raises 
questions concerning their structure, geometry, and dynamics. 


Prospective NASA TN in preparation. 



There are several other questions of intrinsic scientific interest. For 
instance, since Saturn is so similar to Jupiter in other respects, it may be 
expected to have a large magnetic field as well. This raises questions con- 
cerning the interaction between the rings and this possible magnetic field. 
Also, the surface roughness of the particles in the rings may give an indica- 
tion of the average micrometeoroid flux in the vicinity of Saturn. Any incon- 
sistency between this, flux and the present one measured in situ on the same 
flight would suggest that the micrometeoroid flux has varied significantly in 
the past. 

The properties of the individual ring particles are of cosmological inter- 
est. As several theoretical studies have pointed out (refs. 1-3), the rings 
occupy a region of space where, due to tidal forces, the presatellite material 
cloud would not be able to condense into a satellite. If this interpretation 
of the origin of the rings is correct, the rings may represent one of the very 
few regions of the solar system where the remains of the preplanetary cloud can 
be found and studied. Also, the dynamics of the rings' particles may be very 
similar to those of the preplanetary material cloud. 

Because of the limitations of Earth-based observations, very little can 
be said concerning the properties of individual ring particles; in fact, only 
composition and size have been dealt with to date. Analysis of the spectrum 
of reflected infrared solar radiation allowed Pilcher et at. (ref. 4) to iden- 
tify water ice as a major constituent of the ring particles. Using high- 
resolution spectrophotometry, Kuiper et at. (ref. 5) have identified six 
absorption bands of water ice in the light from the rings. Cook 'et at. (ref. 6) 
suggest that impurities, such as dust, are necessary in order to redden the 
visible spectrum. Lebofsky et at. (ref. 7) also contend that the ring parti- 
cles are not pure water ice and suggest the possibility of frost-covered 
silicates . 


As shown in table 1, recent estimates (refs. 8-12) of the size of the 
individual particles have varied over a range of nine orders of magnitude 
(0.1 ym to 100 m) . The hypothesis of small particle sizes, however, raises 
questions of the stability of the rings. Radzievsky (ref. 13) has considered 
the effect of radiation braking on particles in circular orbits about a planet. 
For a sun-planet distance of R (AU) , he found that the time required for the 
orbit of spherical particles of radius a (cm) and mass m (gm) to decay from 
radius r Q to radius zy is given by 


t = 2.3xl0 6 (R/a) 2 m r ln(r o /rj,) years (1) 

The lifetimes given in table 1 assume an average value of R for Saturn of 
9.54 AU, a particle density of 0.9 gm/cm 3 (ice), r = 1.37xl0 5 km (outer radius 
of ring A), r ^ - 9><10 4 km (inner radius of ring B) , and an additional factor 
of 3 to account for the effects of shadowing by the planet and other ring par- 
ticles (ref. 11) . 

If the particles which comprise Saturn's rings are spherical, the aline- 
ment of their spin angular momenta relative to the ring plane is immaterial to 


2 



a consideration of the optical properties of the rings. In the more likely 
case of nonspherical particles, however, the geometrical alinement of the par- 
ticles must be considered. In this study, we consider the behavior of spin- 
ning nonspherical particles in the presence of a magnetic field. The previous 
treatments of this topic have been associated with attempts to explain star- 
light polarization observations in terms of the alinement of interstellar 
grains. The pioneering work in this field was done by Davis and Greenstein 
(ref. 14), (with an addendum by Davis (ref. 15)). At the same time, Spitzer 
and Tukey (ref. 16) proposed a different alinement mechanism, but it has 
largely been ignored since it is much less efficient in small magnetic fields 
than that proposed by Davis and Greenstein. The Davis and Greenstein analysis 
has been put on a firmer theoretical foundation and has been extended to other 
magnetic alinement mechanisms by Jones and Spitzer (ref. 17) and Henry (ref. 18; 
alinement mechanism only) . Work in this area has been briefly reviewed by 
Dieter and Goss (ref. 19); recent work has included that by Greenburg (ref. 20), 
Purcell (ref. 21), Greenburg and Shah (ref. 22), Shah (ref. 23), and Goldstein 
(ref. ,24). ■ 

These previous theoretical analyses have dealt almost exclusively with 
idealized models based on the assumptions of nearly spherical particles and a 
weak magnetic field. Although it is probably not unreasonable to assume that 
the ring particles are nearly spherical, the magnetic fields that could be 
expected in the rings are four or five orders of magnitude greater than those 
in interstellar space. Therefore, we will investigate the degree of alinement 
of nearly spherical particles at equilibrium in the presence of a strong mag- 
netic field. Since the processes that influence the orientation of a real par- 
ticle in a magnetic field are, in general, stochastic in nature, the analysis 
of this problem is centered around the Fokker-Planck equation. Contributions 
to the diffusion coefficients from each of the major factors influencing par- 
ticle orientation are evaluated and approximate solutions of this equation at 
equilibrium are presented. The implications of these solutions for the Saturn 
ring system are then analyzed. 

In order to simplify the discussion of the optical properties of the rings 
in a subsequent study in this series (Part III), we will confine our current 
analysis to the orientation of spheroidal ring particles. 


GENERAL SPHEROID ALINEMENT 


Consider a typical spheroidal ring particle, whose shape will be charac- 
terized by a symmetry semi-axis, a s , and a transverse semi-axis, ay. The par- 
ticle will be assumed to be spinning with an angular velocity, u> (angular 
momentum, ^) ; . the local magnetic field is denoted by B. Let (x, y , s) and . 
(5, n, C) be orthogonal coordinate systems (fig. 1) defined so that e z = eg 
and e g = e s ( e s = unit vector along the particle symmetry axis) . If we let 

be the fraction of spheroids with angular momentum around the i-axis 
between J y and J ^ + dJ y, we can express the Fokker-Planck equation as 
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( 2 ) 




• _ x, y, s 

5, n, c 


where <AJ^> represents the average change in J £ per unit time, and where 
< (AJ) | > represents the average value of (Ad)! per unit time; these quantities 
are referred to as "diffusion coefficients." If we denote by< < AJ^>^, 
<A and<A<7^>y the value of this diffusion coefficient resulting from col 
lisions with other ring particles, from collisions with "dust" (or gas) 
particles, from magnetic torques, and from fluctuations due to the internal 
temperature of the ring particle, respectively, then from the linearity of 
equation (2), we have 


<A J.> = <A J.> 

■z. ^ r 


<*i>d * - <hJ i>T 


. X, 


y> 

n. 


(3) 


A similar notation has been adopted for the other diffusion coefficient as 
well. 


Gas/Dust Collisions 

2 

Expressions for<A and < (AJ) can be derived in a straightforward 
manner (ref. 24). For simplicity, we will use the terms "dust" and "dust par- 
ticles" to refer to interplanetary dust particles, gas atoms, and charged 
particles. 

The average number of collisions between dust and ring particles that 
affect J ^ is + a^v ^) , where is the number density of dust par- 
ticles; = a^; = a s ; and v ^ is the average dust particle velocity 

along the i-axis. The average number of collisions affecting J ^ and J ^ can be 
expressed similarly. Since a ring particle can be assumed to be much more 
massive than a "dust particle," then in kinetic equilibrium its translational 
velocity will be much less, and is therefore ignored. 

The quantity <A J{>d has two components: one due to the angular momentum 

of the dust at impact and the other due to the momentum imparted to the dust 
particles as they leave. The first component will be zero, since, on the aver 
age, as many dust particles will strike the spheroid on one side of each axis 
as on the other, even if the velocity field is anisotropic. One of the sim- 
plest assumptions for the other component is that, on the average, as the dust 
particles leave the spheroid, they have the same velocity as the surface ele- 
ment of the spheroid from which they left. Other assumptions are also plaus- 
ible and might lead to slightly different numerical results. Under this 
assumption, however, 

AJ” .^> i , • • = w -rrija* 2 (4) 

\ z-^one collision z a ^ 
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where m d is the average mass of a dust particle and where a \ is the root mean 
square distance from the -i-axis, which can be shown to be 


* 


a 





(5) 


4 wa 

a = S(e) 

3tt 


( 6 ) 


where 2?(fc) is the complete elliptical integral of the second kind, e is the 
eccentricity of the ellipse formed by the intersection of the spheroid surface 
and a plane which includes the symmetry axis, and 


w = 


1 prolate 

aja oblate 
t s 


(7) 


If we assume an isotropic Maxwellian distribution for the "dust particles" with 
an equivalent dust temperature T^, we obtain 




16a s a^y ^{2vhn ^ T ^) 
31 


( 8 ) 


< A J t>d =<AJ n >d = ' 


3 ttV 


(9) 


where k is the Boltzmann constant, and J. is the moment of inertia around the 
• . % 

^-axls . 


The dominant contribution to < (AJ ),">7 is the angular momentum transferred 


J H'd 

to the spheroid by the dust particle in the collision. This leads to 


<^\>d = 


Z2a l( a l + 2a l)( 2a s + a t) V d{ 2mt ^ T df n 


< + 2 (a 2 * a^ Z 
t vs i 


if 
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( 10 ) 


<W\> d = <i^\> d 


Z2a t( a s * a l) 2 ( a2 s + 2a l)( 2a s * a t ) ' , d( 2m, d k3T d) 
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Magnetic Torques 


Davis and Greenstein (ref. 14) dealt with the effects of magnetic torques 
on a spinning spheroid. They pointed out that if £he spheroid is slightly 
paramagnetic, the component of the magnetization, M, perpendicular to ? will 
fluctuate as the spheroid spins. Since the particle is nonspherical, the sym- 
metry axis will nutate around 3 as the particle spins; and, if the magnetiza- 
tion is averaged over one nutation of the particle, there is a nonzero resul- 
tant magnetization perpendicular to J and proportional to x"/w, where x" is the 
imaginary part of the magnetic susceptibility, x^ : 

X M = x' + £x" (12) 

This nonzero magnetization results in a dissipative torque, which tends to 
decrease the magnitude of the components of ? perpendicular to As a result, 
tends to become alined in the same direction as The results of the anal- 
ysis of Davis and Greenstein (ref. 14), in the notation we are using here, can 
be expressed as 


<' 4 *V5? ■ 0 


(13) 




- Jri'B 2 sin 6 

t> 


[«*■♦ a t) cos e J + 2a l sin 6 j] (14) 


<W 2 >M - 


-2 VX”B 2 J Z sin 2 6 


2J a^a) 

S £ 


— jj^a 2 + a 2 


M. . I W J V r , — / -A. A 4 V/ T 

s tf J t J 


a 2 + at jcos 0 T + 2at sin 6 


r] (15) 


where 


cos = ej • e B = J/J 

(16) 

cos 0 T = e T • e = J/J 
J J s C 

(17) 


Internal Thermal Fluctuations 

The fluctuations due to the internal temperature of the grains have 
recently been investigated by Goldstein (ref. 24). There are many mechanisms 
which may contribute to the internal temperature, Tj, of the grain, including: 
collisions with gas and dust particles, the dissipative torques involved in 
the paramagnetic relaxation process, bombardment by low-energy cosmic rays, by 
the remnants of the solar wind plasma, and by any trapped radiation, and energy 
absorption from the radiation field (sunshine, planet shine, starshine, and 
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moonshine). Since the radiation field, will be the dominant process for estab- 
lishing Tj, we may treat Tj as an independent parameter. 

If Tj is nonzero, the magnetization, M, will fluctuate (ref. 17) in a 
random manner. The fluctuating magnetization wi^l in turn produce a randomly 
varying torque perpendicular to B (torque = M x B) , which will in turn produce 
variations in If we represent the random fluctuations in the magnetization 
by AM, then the symmetry of the problem implies 

<(A M) > =0 i = x, y, z (18) 

Is 

In addition, we can assume that the fluctuations of the orthogonal components 
are uncorrelated: 

< (AM) . (AM) .> = 0 . i / j; i,j = x, y, z (19) 

Since changes in J are directly proportional to the applied torques, equa- 
tion (18) implies 

<(A J j)> T = 0 i = x, y, z (20) 

while the fact that the torques are perpendicular to B implies 

<(A J) 2 z > t = 0 (21) 

It can be shown (ref. 24) that the other two diffusion coefficients can be 
expressed as ^ v 

CCA J)l> T = <^1> T = Stt alakljBW'/Su (22) 


Ring Particle Collisions 


If we let n ( [ i ] , [4. 2 ],*) dt be the probability that two ring particles 
with distribution function parameter sets [<j>i] and [<j> 2 ] collide in the time 
interval from t to t + dt, we can express the effects of such collisions on 
the angular momentum distribution in terms of n. To a first approximation, 
however, we will follow the lead of most of the ring stability analyses 


(refs • 

in i? 

and assume 

that 


and 

<(A^)> r « 

<(A 


<(Ac np r « 

<M\> d 


% = 


y, 
K, n. 


(23) 


Fokker-Planck Equation Solution 

In principle, at least, we have now specified enough information to solve 

the Fokker-Planck equation, equation (2), for f .(J ,0> T ,0 T ) , J = obi, pro. 

J J d 
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Since we are interested in the orientation of the ring particles in equilib- 
rium, we need to find a solution to the homogeneous partial differential equa- 
tion formed by setting the right-hand side of equation (2) equal to zero. 

Such a solution can be approximated in the case of nearly spherical spheroids 
(ref. 24). Using some of Goldstein's (ref. 24) notation and modifying his 
results to correspond to our current problem, we obtain the following solution 
for nearly spherical prolate particles: 


•^pro^’W 




= exp < - I cos 


c£ 2 


sin z g T + 


(1 - Q 2 ) 

-22 
2q o i 


cos z 



(24) 


where d\ and e 2 are constants defined by 


and 


Cl 


1 + bTj/T d 

1 + b 


q 2 ( 3 + AbTjVTj) + 1 
q 2 {Z + 4b) + 1 


(25) 


(26) 


and where b is a measure of the relative importance of magnetic torques and 
gas/dust collisions: 

b = (x"/u)VB 2 /2G\> d (2mi d k3' d ) 1/2 (27) 

and b » 1 indicates that magnetic alinement dominates over the effects of the 
interplanetary medium. In a similar manner, bTj/T j is a measure of the rela- 
tive importance of magnetic torques and fluctuations due to the internal tem- 
perature of the grain. In equation (27), G is a constant related to the size 
and shape of the spheroid, and is defined as follows: 

Gpro = a s a \[2qe 2 + qe + (4e 2 - l)sin -1 (e)]/2e 3 (28) 

= a^[2q l +e 3 - q 2 e + (4q 2 e 2 + l)sin“ 1 {qe) ]/2e 3 (29) 

where 

? " a i/ a e 

q q 

It should be noted that Gp ro /a s a| and ff 0 bi/aja£ are functions of q only; these 
two functions are illustrated in figure 2. The approximate solution for the 
case of nearly spherical oblate spheroids is 

4bl (J ’ W ~~ ex P r|" (1 + q2) COs2 + ° l * Sin2 r.n 

q 2 - 1 . , i J 2 ) 3 

°i J l 8 
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The degree to which the particles are alined with the magnetic field can 
be expressed in terms of the average values of (e s • eg ) 2 and (ej • eg) 2 , 
which indicate the degree to which the symmetry axis and the angular momentum, 
respectively, are alined. Therefore, we are interested in finding 


c 


< cos 2 J | d ( cos e T )cos 2 g T j d(cos 6 T ) | dJ 


r j d( cos Qj)J 


8 q 2 kT I 
o g n 


and 


<cos 2 0 ^ = 


i = Obi, pro 


• = — I d(cos g T ) <£( cos 6 T ) I dJ 

1 1 J.i J J 0 


(32) 


( 


1 - cos 2 g T - cos 2 0 


9 n 


3 cos 2 cos 2 i = obi, pro 


(33) 


with f-i normalized so that the integral of over all of phase space is unity. 
Using equations (24) and (31) in equations (32) and (33) , we can obtain for 
Tj < T £ (see next section for justification) : 


<cos 2 


0 

J pro 


1 - (l/il/ i c? 4 )tan“ 1 (cyc 5 ) 


(34) 


< cos 2 


e> 

pro 


2 2 

q 

2 2 2 
e c 2 °b N i 




(35) 


< cos2 6 />obl = 44 [ X - iff sin * 1 ( 03 )] 
C 8^1 


<cos 2 0> 


obi 


N 2 q i{ q2 + q i 0 l) “ ^ q \ G l^ 1 
- _____ 

N z qie q c 3 


(36) 


(37) 


where 


q\=q 2+ 1 (38) 

and where the constants c 3 , . . . , c 9 have been introduced to simplify the 
above expressions : 
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( 39 ) 


,2 = 

'3 


,2 = 


1 - oi 


(40) 


qI. ^2 

o - 2 e 4 S 
~ - ■ 

5 _ 2 2 


2e^ 


(41) 


,2 = 


. 2 2 
2c x q 


(e 2 + 2g 2 c 2 ^ 


(42) 


c 2 = 


e 2 - c 2 
e 6 C 5 


(43) 


e 2 o 2 


(44) 


2 2 

<? 2 = s_£_ 

9 " 2 2 

«1 C 3 


(45) 


N\ and N 2 are the normalization constants for f^ T0 and respectively: 




\qS{ 2 ) 


Jin 


<? + c 

6 5 


(46) 


ec n 


*2 E 


tan -1 (Cg/^) 




8^1 


(47) 


We can, in a natural manner, define angles that are characteristic of the 
degree of alinement: 



cos -1 [(< cos 2 Bj^) 1 ^ 2 ] 


£ 5 pro, obi 


(48) 


10 



(49) 


0 obl E c°s' 1 [(<cos 2 0> obl ) l/2 ] 

®pr° ^ 11/2 ‘ Q pro 5 si ^" 1 [«cos 2 Q> pro ) l/2 3 

Figure 3(a) shows contours of constant in (<?> &/<?) -space for the range 

q e [1.0, 10.0]. Similar results for 0 o bi are shown in figure 3(b), while 
contours of constant $j } -p ro and 0p ro in ( q , bq ) -space are shown for the range 
q e [0.1, 1.0] in figure 3(c) and F (d). 

It is clear that, for large b, obl , 3j^p ro , and 0 obl tend to zero 
while ©p ro tends to tt/ 2. This means tftat, if magnetic effects are dominant, 
oblate spheroids will tend to aline themselves with their spin axes and sym- 
metry axes coincident and parallel to S. Prolate spheroids, on the other 
hand, will tend to aline themselves with their spin axes parallel to B and 
their symmetry axes perpendicular to B. Thus , both types of spheroids tend 
to rotate around the shortest body axis and to aline this axis with a strong 
magnetic field. 


ALINEMENT IN SATURN'S RINGS 

In order to apply these results to the particles in Saturn's rings, we 
must estimate the pertinent particle, planetary, and interplanetary parameters, 
namely, q, a g , T Jt X"/w, B , v^, m^, and 2^. 


Cis-Saturn Parameter Values 

q.— In our discussion of ring particle alinement we will be principally 
concerned with nearly spherical particles {q » 1). As a consequence, although 
we will parameterize our results with respect to q, we will limit q to the 
range [0.1, 10.0]. 

a s .— The range of particle sizes, characterized by a s , is constrained by 
the current estimates of ring particle sizes (see table 1) : 

10" 8 m < a < 10 3 m (51) 

s 

Tj.— The best current estimates of the temperature of the ring particles 
are based on radiometric and radio astronomical measurements . Two temperature 
levels appear to fit the observational data: Tj * 30° K and Tj * 65° K 

(ref. 11) . Since the effectiveness of the magnetic alinement mechanisms is 
inversely proportional to the internal temperature of the particles, we will 
choose the latter level (65° K) as representative of "worst case" conditions. 

X"/a).— The magnetic susceptibility of the particles depends on their 
internal composition, structure, and thermodynamic state. The value to be 
used for (X"/w) also depends on the magnetic relaxation mechanism being 
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considered. Since it does not appear feasible to deal quantitatively with 
diamagnetic effects, as, for instance, in the case of graphite flakes imbedded 
in ice (ref. 29), we limit our considerations to paramagnetic and ferromagnetic 
effects. The derivations above are, in part, based on the assumption that the 
quantity X"/w is approximately constant. This assumption must be justified 
and the constant evaluated; this is done briefly for paramagnetic and ferro- 
magnetic interactions. More detailed discussions are given by Davis and 
Greenstein (ref. 14) and Jones and Spitzer (ref. 17). 

The dissipation of energy by paramagnetic forces is usually described in 
terms of a combination of spin-lattice and spin- spin interactions. However, 
if the paramagnetic ions are homogeneously distributed throughout the particle, 
there is no net magnetic moment in a coordinate system rotating with the par- 
ticle; hence, spin- lattice interactions do not contribute to paramagnetic 
relaxation. The standard method of calculating the effects of spin-spin 
interactions is that of Van Vleck (ref. 30), which results in very tedious 
calculations. In order to avoid these calculations, Locher and Gorter (ref. 31) 
proposed an interpolating lineshape function which has been used successfully 
and predicts 

X"(<jj) = X(o)@ l/2 [a/^ exp (a 2 ) erfc (a)] -1 

exp(-a 2 T 2 co 2 ) (52) 

1 + ai T 

where a is a function of<u) 4 >/< w 2 > 2 and is a measure of the effectiveness 
of such other effects as hyperfine, exchange, and crystal field interactions, 
and x is the time constant for spin-spin interactions. For any reasonable 
values of a (>10~ 2 ) it is easy to show that the expression for X"(w) reduces 
to a Gaussian (within a factor of 2) : 

X f » = X(o)(f) 1/2 - exp (-co 2 /2 <u) 2 >) (53) 

V2/ < a ) 2 > 1/2 

This, of course, leads to the conclusion that X"(to)/w is approximately constant 
for small oj. Jones and Spitzer (ref. 17) have shown that 10 -12 /2 j should 
represent a realistic lower limit for X"/w. They have also pointed out that 
if there are abundant concentrations of hydrogenic atoms, nuclear paramagnetic 
effects obviate the presence of paramagnetic ions, and the lower limit of 
X"/oj is unchanged. Thus, we have 


( x "/w)para * 10 ~ l2 fT J . = 1.5xl0" 14 sec/rad (54) 

nearly independent of the chemical composition of the ring particles. 

Jones and Spitzer (ref. 17) have pointed out that if, instead of the 
uniform, diffuse distribution of paramagnetic impurities assumed above, one 
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considers the equally likely possibility that these impurities may aggregate 
into particulates of, say Fe 3 04 or y-Fe 2 0 3 , then the value of (X"/w)p ara 
might be enhanced by six or seven orders of magnitude. Jones and Spitzer refer 
to this mechanism as "super-paramagnetism:" 

( x '» s _para ' 10 " 8sec/rad (55) 

A consideration of ferromagnetic effects is complicated somewhat by the 
fact that there is a critical size (-0.01-1.0 pm) below which the ferromag- 
netic grains in the ring particles would consist of essentially a single 
domain, and above which they must be treated as being multidomained. Theo- 
retical treatments of multidomain particles are rather sparse, but it is gen- 
erally concluded that the low-frequency behavior of their permeability is due 
to reversible domain-wall motion (refs. 32-34). By appealing to the available 
experimental data, Jones and Spitzer (ref. 17) conclude that X"/u) is constant 
in the limit of small frequencies. Excluding more exotic compounds such as 
Ni 0- 4 Zno.6 Mn o.02 Fei _9 0 4 , the range of values of X"/w for multidomain fer- 
romagnetic particles at room temperatures is approximately (ref. 35) 

10" 11 < ( x "/w) ferr 0 ,md < 10 " ? sec/rad (56) 

Since the temperature dependence of this value is approximately like 
exp(AE/kTj) (0.1 eV < A E < 1.0 eV, ref. 36), we would expect X"/m to be much 
larger than indicated in equation (56) . Experimental data at low temperatures 
is virtually nonexistent, however, and we therefore have used the vastly under- 
estimated range indicated in equation (56) , which results in an underestima- 
tion of the degree of particle alinement. The analysis of the permeability of 
single-domain ferromagnetic particles in the low-frequency limit is reviewed 
by Jones and Spitzer (ref. 17), who conclude that for interstellar grains 

2xl0 -13 < (X"/u)^ erro s( j < 4><10” 7 sec/rad (57) 

"with the lower limit extremely unlikely." 

Clearly, the wide range of possible particle compositions and the several 
different possible relaxation mechanisms precludes a definitive estimate of 
X"/w; it is possible, however, to assert that 1.5x10“ 14 sec/rad is definitely 
a lower bound and that the actual value is probably several orders of mangi- 
tude larger. 

Haffner (refs. 37, 38) has investigated the possibility of a 
Saturnian magnetic field. On the basis of the limited radioastronomical data 
available for Saturn and the assumption that the rings will effectively sweep 
trapped charged particles out of the inner magnetic field, he estimates the 
magnetic field strength at the equator to be 1 to 10 gauss. If we make the 
reasonable assumption that the magnetic field near the equatorial plane is 
essentially dipolar out to at least 2.5 Bp,, then the direction of the magnetic 
field in the rings will be normal to the ring plane and the magnitude of the 
field will vary as pp, -3 , where pp, is the radial distance from the center of 
the planet. Table 2 shows the range of magnetic field strengths to be 
expected at the boundaries of the major ring elements. 
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v d> m d* T d'~ T ^ e characteristics of the interplanetary dust/gas environ- 
ment have been estimated by Cook and Franklin (ref. 27) in connection with a 
study to determine the effect of this environment on the stability of the 
rings. They estimate that the mass density of dust particles in cissaturnian 
interplanetary space is about 1.2xl0 -23 gm/cm 3 . If we assume that the density 
of the individual dust particles is the same as that of meteoroids, 0.18 gm/cm 3 
(ref. 39), then, assuming a dust particle radius of 1.0 ym, we have 

v d Z 1.6X10" 12 cm -3 (58) 


and 


m d = 7.5xl0 -12 gm (59) 

Cook and Franklin also estimate that the average dust particle velocity is 
about 34 km/sec, which leads to an effective kinetic temperature: 

T d z 2.5xl0 17 °K (60) 

Degree of Ring Particle Alinement 

With these estimates of the pertinent particle, planetary, and inter- 
planetary parameters, we can now estimate the degree to which the particles in 
Saturn's rings could be expected to be alined. If we assume specific values 
for 5 f>, equator and X"/oj then we can express this degree of alinement in terms 
of the maximum size particle that would be alined to within a specific angle 
as a function of position in the rings for each of the four types of magnetic 
relaxation methods. Figure 4 (a) -(d) shows such a display for the two aline- 
ment angles for each of the two types of spheroids. Here we have chosen 
equator = ^ gauss as typical, and for each of the four relaxation methods 
we nave used a value of X"/ui that is two orders of magnitude above the abso- 
lute lower bound. All these choices appear to be conservative, and more real- 
istic estimates of a Sjmax are probably several orders of magnitude greater. 

The data in figure 4(a), for instance, indicate that, with these para- 
metric estimates, the spin axes of all prolate spheroids in the B ring smaller 
than 3xl0 -6 cm would be alined to within 1° of B even for the least efficient 
relaxation method (paramagnetic relaxation); if multidomain ferromagnetic par- 
ticulates occur in the particles, however, their spin axes would be alined to 
within 1° of B for a s < 2.5 cm. The effect of uncertainty in the correct 
value for X"/to, is indicated by the use of the upper bound for this parameter, 
which changes these two values to a s < 3xl0 -2 cm and a s <_2.5 m. Similar con- 
clusions can be drawn from figure 4(b), (c) , and (d) for 0p ro , $J 3 obi ’ anc * 

0 bi . Fully alined particles would appear as shown in figure 5. 3 
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CONCLUSIONS 


If the particles in the rings of Saturn are spherical, the orientation of 
their spin angular momenta is immaterial to their optical properties. Per- 
fectly spherical particles are, however, quite unlikely, and the orientation 
of nonspherical particles can be anticipated to have a significant effect on 
the optical properties of the rings. As a basis for subsequent investigation 
of these effects (Part III), we have determined the "preferred" orientation of 
spheroidal particles in a magnetic field and investigated the degree to which 
the ring particles will assume this orientation. 

Since the alinement mechanisms, as well as the misalinement mechanisms, 
are stochastic in nature, we have used the Fokker-Planck formulism to deter- 
mine the alinement distribution function. The four major factors contributing 
to the particle orientation — the interplanetary dust/gas, the Saturnian mag- 
netic field, the internal ring particle temperature, and the collisions 
between ring particles — have been analyzed in turn, and the contribution of 
each to the Fokker-Planck diffusion coefficients has been determined. The 
solution of the resultant Fokker-Planck equation for nearly spherical par- 
ticles has been presented, and the particle alinement implied by this solu- 
tion (fig. 3) indicates that, at equilibrium, in the presence of a suffi- 
ciently large magnetic field any spheroidal particle will be spinning around 
its shortest body axis, with its spin axis oriented parallel to the magnetic 
field. 

From estimates established, for the pertinent physical parameters in the 
vicinity of the rings, implications of the previous analysis for the ring par- 
ticles were investigated (fig. 4). The significance of the resulting data is 
somewhat obscured by the uncertainty in the value of the imaginary part of 
the magnetic susceptibility, x"/u). Nevertheless, it is clear that alinement 
of the ring particles is quite likely, especially for small (<10 cm) particles 
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TABLE 1. - ESTIMATED SIZE OF SATURN RING PARTICLES 


Source 

Characteristic radius. 

Radiation braking lifetime 

meter 

for A § B rings 

(after ref. 13) 



Years 

Solar system 



lifetimes 

Schoenberg 
(ref. 8) 

- 3.6 x 10" 6 

* 3.5 x 10 5 

> 7.8 x 10“ 5 

Franklin and 
Cook (ref. 9) 
("Model II") 

3.1 x 10" 4 

3 x 10 7 

6.7 x 10‘ 3 

Franklin and 
Columbo 
(ref. 10) 

1 x 10 2 

1 x 10 1 3 

2.2 x 10 3 

Bobrov 

1 x 10 -6 - 7.5 

1 x 10 5 - 7.3 x 10 11 

2.2 x 10' 5 - 1.6 x 10 2 

(ref. 11) 

4 . 3 x 1 0“ 2 - 2 . 8 x 10 2 

4 . 2 x 10 9 - 2.7 x 10 13 

9.3 x 10' 1 - 6 x 10 3 

Price 
(ref. 12) 

1 x 10" 7 

1 x 10 4 

2.2 x 10" 6 

\ 

Pollack 

(unpublished 

estimate) 

6 x 10~ 5 - 1.0 

6 x 10 & - 1 xlO 11 

1.3 x 10~ 3 - 22 

Minimum radius 
for lifetime 
comparable to 
solar system 
lifetime 

4.6 x 10" 2 

4.5 X 10 9 

.... . 

1.0 
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TABLE 2. - DIMENSIONS OF SATURN'S RINGS (ref. 42) 



Ring boundary 

Ring boundary 

Magnetic field 


radius ( 10 3 km) 

radius (/?,) 

strength (est.. 




ref. 38) (gauss) 


270? (ref. 41) 

4.5? 

0.010 - 0.10 

Ring D 

137 

2.28 

0.084 - 0.84 

Ring A 

120 

2.00 

0.125 - 1.25 

Cassini division 

116 

1.93 

0.139 - 1.39 

Ring B 

90 

1.50 

0.296 - 2.96 

Gap 

89 

1.48 

0.308 - 3.08 

Ring C 

72 

1.20 

0.579 - 5.79 

Saturn 

60 

1.00 

1.00 -10.0 
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(c) Prolate spheroids, spin axis. (d) Prolate spheroids, symmetry axis. 

Figure 3.— Alinement of the spin and symmetry axes of spheroids with the mag- 
netic field. The quantity, q, is the ratio of the transverse and symmetry 
semi-axes (<7 = a^/a s ), and b is a measure of the relative effectiveness of 
magnetic and dust/gas interactions; (b » 1 implies that magnetic effects 
dominate) . 
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RELAXATION METHOD 



P'\ 

(a) Prolate spheroids, spin axis. 


RELAXATION METHOD 



. P' R ll 


(b) Prolate spheroids,, symmetry axis. 

Figure 4.— Degree of alinement of the spin and symmetry axes of spheroids with 
the Saturnian magnetic, field.- The abscissa is the radial distance from the 
center of the planet.: in units of Saturn radii (f?^) . The ordinate scales 
represent the maximum "size ring particles (in terms of a g ] , which will give 
the indicated degree of alinement for each of. the four magnetic relaxation 
mechanisms indicated. . 
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RELAXATION METHOD 



obi = 10° 


0j,obl"° 

^j,obr OJ 


(c) Oblate spheroids, spin axis. 


RELAXATION METHOD 



(d) Oblate spheroids, symmetry axis. 
Figure 4.— Concluded. 
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(b) Oblate. 


Figure 5.— Fully alined spheroids. 
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